A silicon-chip source of bright photon-pair comb 
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Abstract 

Integrated quantum photonics relies critically on the purity, scalability, integrability, and flexibility of a 
photon source to support diverse quantum functionalities on a single chip. Up to date, it remains an open 
challenge to realize an efficient monolithic photon-pair source for on-chip application. Here we report 
a device on the silicon-on-insulator platform that utilizes dramatic cavity enhanced four-wave mixing in 
a high-Q silicon microdisk resonator. The device is able to produce high-purity photon pairs in a comb 
fashion, with an unprecedented spectral brightness of 6.24 x 10^ pair/s/mW^/GHz and photon-pair cor- 
relation with a coincidence-to-accidental ratio of 1386 ±278 while pumped with a continuous-wave laser. 
The superior performance, together with the structural compactness and CMOS compatibility, opens up a 
great avenue towards quantum silicon photonics with unprecedented capability of multi-channel parallel 
information processing for both integrated quantum computing and long-haul quantum communication. 
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Photon-based quantum information technology has found broad applications ranging from 
quantum communication \ quantum computing, to quantum metrology^. Recent advances in 
integrated quantum photonics"^"^ show great promise for chip-scale quantum information process- 
ing with tremendous complexity. A bright, single-mode, high-purity, integrated source of single 
photons and/or entangled photon pairs is essential for all these applications'^^^, particularly for in- 
tegrated quantum photonic circuits which rely critically on the purity, scalability, integrability, and 
flexibility of the photon source to support diverse quantum functionalities on a single chipP^^. 
However, although a variety of photon sources have been developed in the past few 

decadePS]^ it 

remains an open challenge to realize an efficient monolithic source for on-chip application, lack of 
which current quantum photonic devices have to rely on external sources for proper operatioiP^. 
Here we propose and demonstrate an ultra-bright high-purity chip-scale photon-pair source on 
the silicon-on-insulator (SOI) platform. By taking advantage of the dramatic cavity enhanced 
four-wave mixing in a high-quality silicon microdisk resonator, we are able to achieve a spec- 
tral brightness of 6.24 x 10^ pair/ s/mW^ /GHz orders of magnitude larger than other photon-pair 
source^^^J^^^, and an unprecedented photon-pair correlation with a coincidence-to-accidental 
ratio of 1386 ±278 while pumping with a continuous-wave laser. In particular, the unique de- 
vice characteristics enable photon pair generation in multiple frequency combs, thus significantly 
extending the wavelength management capability for integrated quantum photonics. The unprece- 
dented device performance together with its CMOS compatibility now opens the door towards 
quantum silicon photonics with great potential for multi-channel parallel operation of novel quan- 
tum functionalities on chip. 

To date, nearly all photon-pair sources are based upon spontaneous parametric down conver- 
sion (SPDC) or four-wave mixing (FWM) in nonlinear optical crystals/waveguideJ^'^. Bulk 
crystals emit photon pairs into a multimode spatial profile, resulting in a fairly low photon gen- 
eration/collection efficienc}^ 1 | 1 1 1 14 | 1 5 J Significant efforts have been devoted in recent years to 
developing waveguide sources for single-mode emissiorPElHIT] However, the produced photon 
pairs generally exhibit a non-factorable spectrum which degrades considerably the quantum-state 
purity^. FWM in silica optical fibers recently appears as a promising approach with a great flexibil- 
ity of engineering photon spectrunP^^lZl^ which, unfortunately, suffers seriously from broadband 
Raman scattering of 

silicdMMI. 

Silicon exhibits a strong Kerr nonlinearity for nonlinear optical interaction and a large refractive 
index enabling tight mode confinement, which has been explored intensively recently for a vari- 
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FIG. 1 : a. Schematic of generating a photon-pair comb from a siHcon microdisk resonator sitting on a silica pedestal. 
All photon modes couple to the same propagation mode inside the delivery coupling waveguide, which is a single- 
mode tapered silica fiber, b. Group-velocity dispersion of vertically polarized transverse-magnetic-like (TM-like) 
cavity modes with four different radial orders, termed as Pq, Pi, P2, and P3, respectively, for a silicon microdisk with 
a thickness of 260 nm and a radius of 4.9 /zm, simulated by the finite element method (Appendix [A|. The insets show 
the simulated optical field profiles. 

ety of applicationd^^^^. In particular, single-crystalline silicon has a clean phonon spectrum with 
Brillouin- zone-center phonons of a well-defined frequency of 15.6 THz and a narrow linewidth of 
~105 GHz, which eliminates the deleterious broadband Raman noise^^mil These superior fea- 
tures together with mature nanofabrication technology make SOI an ideal platform for integrated 
quantum photonic applicatioiP^W35|. 

The device structure we employ for high-quality photon-pair generation is a compact silicon 
microdisk resonator sitting on a silica pedestal (Fig. [T^, with the scanning-electron microscopic 
image shown in Fig. |2]). An excellent feature of a microdisk is that its group- velocity disper- 
sion is dominantly determined by the disk thickness, allowing flexible dispersion engineering for 
the FWM process. As shown in Fig. [I]?, with a disk thickness of 260 nm, the zero-dispersion 
wavelength (ZDWL) for the vertically polarized transverse-magnetic-like (TM-like) mode can be 
tailored to the useful telecommunication band around 1.5 jUm. Detailed analysis shows that the 
dispersion characteristics support phase matching over a broad spectrum around 100 nm, thus en- 
abling simultaneous generation of photon pairs in a comb fashion over all phase-matched cavity 
modes. 
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In particular, the microdisk resonator uniquely supports multiple mode families with different 
radial orders (Fig. [IJ5, inset). As they experience different mode confinement introduced by the 
whispering-gallery geometry, these mode families exhibit ZDWLs shifted by about 30-40 nm from 
each other (Fig.[lj5). Consequently, a variety of desired photon-pair combs can be produced inside 
a single device simply by pumping at different mode families. Moreover, the complete air cladding 
ensures a clean interaction of the optical field with the silicon core and thus eliminates any potential 
contamination from the cladding material, in contrast to waveguides/microrings^^^^^^'^^Hlll Y^,hej-e 
the optical interaction with the buried oxide layer is likely to produce Raman noise photons. 

Most importantly, the high-Q cavity drastically modifies the density of states of vacuum, re- 
sulting in significant Purcell enhancement on cavity modes which not only dramatically increases 
the generation efficiency of photon pairs but also ensures them to be created at discrete single 
frequencies with a high state purit}^. With Np pump photons inside the cavity, the probability of 
emitting a pair of signal and idler photons at time ts and ti into the coupling waveguide (Fig. [T^) is 
given by (Appendix |B]) 

p,(.„.,0 = ^W)Vr..k'-'i (1) 

with Ftj = Fts (or Fti) on the exponent when ts > ti (or t^ < ti), where Fej and Ftj (j = s,i) are the 
external coupling rate and photon decay rate of the loaded cavity, for the signal and idler photons, 
respectively, and F = (Fts -|- Fti) /2 represents the average, g = describes the vacuum 

coupling rate of the FWM process, where n2, V, and rj are the Kerr nonlinear coefficient, the 
effective mode volume, and the spatial mode overlap among the interacting photons, respectively. 

As indicated by Eq. ([!]), the signal and idler photons are created in a time scale of the cavity 
photon lifetime within which the photon pair remains highly correlated. The longer the photon 
lifetime, the larger the pair emission probability. The resulting photon-pair emission rate is given 
by 

= T^-j^ F ■ 

i tsi ti t 

FWM creates photon pairs inside the cavity at a rate of 2{gNp)^ /T, which are delivered to the 
coupling waveguide with a pair extraction efficiency of FesFei/ (FtsFti). A high cavity Q, a strong 
Kerr nonlinearity, and a small effective mode volume would result in a large photon-pair emission 
rate. For example, simulations by the finite-element method show that a silicon microdisk with 
a radius of 5 jum exhibits a small effective mode volume of ~ 7 jUm-^ for a Pq mode, resulting 
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FIG. 2: Schematic of the experimental setup for optical spectroscopy and coincidence photon counting. The pump 
laser is coupled into the device via a single-mode tapered silica fiber. Before coupling into the device, the pump laser 
passes through a bandpass filter and a coarse wavelength-division multiplexer (CWDM MUX) to cut the laser noises. 
The CWDM MUX has a 3-dB bandwidth of 17 nm for each of its transmission bands whose center wavelengths are 
separated by 20 nm apart with a band isolation of > 120 dB. The bandpass filter is identical to one band of the CWDM 
MUX. The produced photon-pair comb is separated into individual photon modes by the CWDM DEMUX which is 
identical to the CWDM MUX used at the input end. The photoluminescence spectrum of the photon-pair comb is 
recorded at each transmission port of the CWDM DEMUX for easy suppression of the pump wave. For coincidence 
counting, the photon pairs are recorded by gated InGaAs single photon detectors (SPDs). Two tunable bandpass 
filters with a 3-dB bandwidth of 1.2 nm are used to cut the Raman noises produced by the delivery silica fiber The 
pump power is controlled by variable optical attenuators (VOAs) and its wavelength is calibrated by a Mach-Zehnder 
interferometer (MZI). The inset shows the scanning electron microscopic image of the device. 

in a FWM coupling rate of ^ ~ 2;r x 460 Hz. Therefore, even a small pump power of 10 /iW 
launched into a critically coupled cavity with a Q factor of 10^ is able to emit photon pairs at a 
rate ~ 6 x 10^ pair/s, clearly showing the extremely high efficiency of the proposed device. 

The fabricated device (see Method) was tested using the experimental setup shown in Fig. [2} 
For easy separation of individual photon modes within the generated photon comb, the device was 
designed with a free spectral range (FSR) close to 20 nm such that the cavity resonances fit directly 
into the transmission bands of a standard coarse wavelength-division multiplexer (CWDM). The 
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FIG. 3: a, The transmission spectrum of the passive cavity scanned by two tunable lasers operating at different 
spectral ranges (indicated as blue and green). The two sets of arrows indicate cavity modes with constant FSRs within 
the Pi and P2 mode families, respectively. The insets show detailed transmission spectra of a P2 mode at 1497.2 nm 
and a Pi mode at 1564.9 nm, with theoretical fitting shown in red. b,c. The photoluminescence spectra of the P2 
and Pi photon-pair comb, respectively, with the pump nearly critically coupled to the cavity. The colors indicate the 
individual spectra recorded at different transmission ports of the CWDM DEMUX. The four mode pairs in the two 
combs are termed as Pj, P2 , P^, and P'\ respectively, as indicated on the figures. The spectrum of the Pi comb extends 
to the photon mode at 1584.1 nm which is beyond our laser scanning range and is thus not seen in a. 

pump wave and the created photon pairs are delivered into and out of the microdisk through a 
tapered fiber which provides flexible control of light coupling. 

The laser-scanned transmission spectrum of the passive cavity (Fig. |3^) shows four clear high- 
Q mode families. A constant FSR of 2.4225 THz (corresponding to ~18 nm around 1.5 jim) 
is observed to maintain among five cavity modes for the P2 mode family, covering the spectrum 
from 1478 to 1555 nm. Similar phenomenon is observed for the Pi mode family with a constant 
FSR of 2.3158 THz, over the spectral region beyond 1510 nm. A detailed scan of the P2 mode at 
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1497.2 nm and the Pi mode at 1564.9 nm (Fig.js^, insets) shows intrinsic Q factors of 3.47 x 10^ 
and 4.94 x 10^, respectively. Cavity modes within a same mode family exhibit similar optical 
Q since, with a same mode profile, they are of similar sensitivity to scattering loss induced by 
fabrication imperfection. These two mode families are used below for photon-pair generation. 
The Po and P3 mode families are not used in experiment since their ZDWLs (~1572 and 1466 nm, 
Fig. [l]?) are outside the scanning ranges of our lasers, although the Pq mode family exhibits a 
higher intrinsic Q factor of about 1 x 10^. 

The constant mode spacing together with a high optical Q ensures efficient generation of a 
comb of photon pairs. As shown by the photoluminescence (PL) spectra in Fig. [3J5, pumping at 
the P2 mode located at 1515.6 nm is able to produce a comb of four clean photon modes within 
the P2 mode family (Fig. [3J5). Moreover, switching to pump at the Pi mode at 1546.3 nm is 
able to produce another set of clean photon comb within the Pi mode family (Fig. [3]:). For the 
simplicity of notation, we term the photon mode pairs in the two combs as F\, Pj^, and P2, 
respectively (Fig. |3j5,c). The spectrum of each photon mode is so sharp that it is beyond the 
resolution of our spectrometer (~ 0.135 nm), implying the high coherence of generated photons. 
The amplitude difference within a photon comb is primarily due to different external couplings 
of cavity modes to the tapered fiber. When the pump mode is critically coupled to the cavity, 
the signals at shorter wavelengths are under coupled while the idlers at longer wavelengths are 
over coupled, resulting in a higher photon extraction efficiency for the idler within each mode pair 
(Fig. [3|5,c). It also leads to a smaller overall loaded Q for the Pj (Pi^) mode pair compared with 
P2 (Pj), thus resulting in a lower pair production rate. Note that the PL spectra also show some 
other modes with small amplitudes whose physical origin is not clear at this moment. Potential 
origins are multiphonon scattering in siliconP'or light scattering from certain defects on the device 
surface, which require further investigation. As these tiny modes are well separated from the 
FWM-created photon combs, they have negligible effect on the comb performance. 

To characterize the performance of the photon-pair combs, we carried out coincidence photon 
counting for each mode pair when the device is pumped with a continuous-wave (CW) laser. As 
shown by Fig. |4^-h, the coincidence counts for all correlated mode pairs depend quadratically on 
the pump power, indicating they are all produced through the FWM process. The experimentally 
recorded pair fluxes (blue dots) are close to the theoretical predictions (red curves) (Appendix |B]), 
with a small discrepancy primarily coming from the uncertainty of Kerr nonlinear coefficient as 
well as measured optical Q. No significant saturation is observed, thus indicating the negligible im- 
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FIG. 4: Recorded photon emission fluxes, pair emission fluxes, and CARs for the photon-pair combs. Each column 
of figures corresponds to one mode pair as indicated on the top. The figures within each row are plotted with identical 
vertical axes for each comparison among different mode pairs, a-d. Recorded emission fluxes of individual photon 
modes within each mode pair, with detector dark counts subtracted. The signal (idler) at the shorter (longer) wave- 
length is shown as blue (red). Solid curves show the theoretical prediction of pure FWM-created photon fluxes, e-h, 
Recorded emission fluxes of correlated photon pairs. In each figure, the left and right vertical axes show the emission 
flux and the corresponding true coincidence count per gate, respectively, indicated by brown and green colors. The 
solid curves show the theoretical prediction, i-1. Recorded CAR for each mode pair No any accidental coincidence 
is subtracted. The red solid curves show the theoretical prediction. The black dashed curves are calculated from 
independently recorded true coincidence counts and photon fluxes of individual modes (Appendix [C|. All data were 
recorded with a detector gating width of 10 ns, a quantum efficiency of 15%, a clock frequency of 250 kHz, and a 
dead time equal to the clock period. 



pact of free-carrier absorption and two-photon absorption. In particular, a small input pump power 
of 79 jUW is able to generate a photon-pair flux of 8.55 x 10^ pair/s at V\ (Fig.|4^), clearly show- 
ing the high efficiency of the device. A metric characterizing the photon-pair generation efficiency 
is the spectral brightness, defined as the pair flux per unit spectral width per unit pump power 
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square (since FWM depends quadratically on the pump power). The spectral width of each photon 
mode is determined by the linewidth of the loaded cavity, which is 1.59 and 2.80 GHz for the 
mode pair at 1497.2 and 1534.3 nm, respectively. It infers that the spectral brightness of the gen- 
erated photon pairs is 6.24 x 10^ pair/ s/mW^ /GHz, which is about 3 orders of magnitude larger 
than any FWM -based photon-pair sources reported up to dateP^^. P\, and Pj^ exhibit similar 
spectral brightness, with magnitudes of 1.78 x 10^, 2.32 x 10^, and 2.57 x 10^ pair/s/mW^/GHz, 
respectively. 

Indeed, the cavity enhancement is so strong inside the device that the demonstrated photon- 
pair source is even brighter than SPDC -based sources, although FWM, a third-order nonlinear 
process, is generally much weaker than SPDC which is a second-order nonlinear process. The 
spectral brightness of SPDC-based sources is generally described by the slope efficiency given as 
the photon-pair flux per unit spectral width per unit powefl^^^. Pj shown in Fig. |4^ has a slope 
efficiency of 4.93 x 10^ pair/s/mW/GHz at a pair flux of 8.55 x 10^ pair/s, which is more than 
one order of magnitude larger than any SPDC-based waveguide source^^lMEol jj jg comparable 
with the state-of-the-art cavity-enhanced SPDC sourced^^^^, while the latter create photon pairs in 
multi-spatiotemporal modes which require stringent filtering. Moreover, the latter are constructed 
with bulk optical components which require sophisticated cavity locking and are challenging for 
chip-scale integration. 

Figure |4] shows that the pair fluxes are smaller than those of corresponding individual photon 
modes. This is because of the non-unity photon extraction efficiency, a common feature of cavity 
quantum electrodynamic systems: Although the signal and idler photons are always created in a 
pair fashion inside the cavity, they transmit out of the cavity independently, only a certain fraction 
of which remains the pair correlation. The performance can be improved by using cavity modes 
with higher intrinsic Q (e.g., the Pq mode family) together with higher external coupling. On the 
other hand, Fig.|4^-d show that the fluxes of individual photon modes deviate slightly away from 
the quadratic power dependence, indicating that they are accompanied with a certain amount of 
noises. Detailed characterizations show that the noise photons dominantly come from the broad- 
band Raman noises generated inside the delivery fibers (Appendix [P]). In experiment, a pair of 
bandpass filters with a 3-dB bandwidth of 1 .2 nm were used in front of the single photon detectors 
(SPDs) (Fig. [2]) to suppress the Raman noises. As each photon mode has a very narrow spectral 
width, the Raman noises can be suppressed further by using a narrower filter. They can also be 
reduced by shortening the lengths of delivery fibers. Such Raman noises would be absent if an 
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on-chip lasei'^Slis used to generate the photon-pair combs. 

An important figure of merit characterizing the quantum correlation between the photon pairs 
is the coincidence-to-accidental ratio which is directly related to the quantum- 

interference visibility of the entangled states constructed from the photon pairs^^. Figure Q-l 
show the recorded CAR (blue dots) for all mode pairs. In general, CAR increases with decreased 
photon-pair flux because of reduced probability of multi-pair creation and saturates at a certain 
level of pair flux when the detector dark counts start to dominate. All mode pairs exhibit very high 
CAR at all power levels. For the P2 mode pair, the device is able to achieve unprecedented CAR 
with a value of 274 ± 10 at a pair flux of 8.55 x 10^ pair/s and a peak value of 1386 ± 278 at a 
pair flux of 2.35 x 10^ pair/s. These values are more than three times larger than other photon- 
pair source^ even without using pulsed pump to suppress noise photon^ ^ l ^^^ l ^^ l ^^ l ^^ l ^ ^ l ^^ l or 
dehcate superconducting SPDs to reduce detector dark countd^ ^^ ^ ^^ l ^^ l Similar CAR is achieved 
for P{, with values of 1097 ± 165 and 244 ± 6 at pair fluxes of 2.43 x 10^ and 6.31 x 10^ pair/s, 
respectively. At the same time, and Pj exhibit as well relatively high CAR, with peak values of 
733 ± 106 and 679 ± 74, respectively. For all the mode pairs, the power-dependent trend of CAR 
follows closely the theory (red curves) (Appendix |C]). The discrepancy primarily comes from the 
uncertainty in characterizing Raman noises generated in delivery fibers. These results clearly show 
the superior quality of generated photon-pair combs. We expect that the peak CAR of the photon- 
pair combs can be improved by at least one order of magnitude in the future, by suppressing the 
Raman noises produced in the delivery fibers and by using better detectors (e.g., superconducting 
SPDs) with lower dark counts. 

The high pair correlation and high coherence of the photon-pair combs allow resolving the 
temporal structure of quantum correlation between the signal and idler photons. This is particularly 
enabled by the tapered-fiber coupling system which provides flexible control of the external light 
coupling and thus the cavity photon lifetime. Figure [5] shows the normalized coincidence spectra 
of P2 with two different external coupling conditions. When the tapered fiber is close to the device, 
the external coupling is strong which results in a short photon lifetime of 65 ps (averaged between 
signal and idler). As a result, the coincidence spectrum (blue curve) shows a full width at half 
maximum (FWHM) of 370 ps, which is 19% larger than the time resolution of the instrument 
response function (IRF, gray curve) of our coincidence counting system. However, moving the 
tapered fiber away from the device increases the photon lifetime considerably to 171 ps, leading 
to a coincidence spectrum (green curve) with a FWHM of 453 ps which is 45% larger than the 



10 




Delay Time (ps) 

FIG. 5: Normalized delay-time-dependent coincidence spectra for Pj. The blue and green curves show the cases 
with different external coupling conditions and thus different averaged photon lifetimes. The gray curve shows the 
instrument response function (IRF) of our coincidence counting system, with a FWHM of 312 ps which primarily 
come from the timing jitters of the two SPDs. The red and purple curves show the theoretical prediction. The 
coincidence spectra were recorded with a detector gating width of 10 ns, a quantum efficiency of 25%, and a clock 
frequency of 125 kHz. The coincidence counter has a time resolution of 4 ps. 

time resolution of IRF. The experimentally observed time correlation is explained well by the 
theory (Appendix |C]). To the best of our knowledge, this is the first time to resolve the temporal 
photon correlation for a chip-scale photon-pair source, although it was observed before in cavity- 
enhanced SPDC source^^^^. Quantum coherence between correlated photon pairs is essential for 
constructing a variety of quantum functionalitie^^^^. The photon lifetime of 171 ps in our device 
corresponds to a propagation length of ~ 1.5 cm, which is about the footprint of a photonic chip. 
Therefore, the produced photon pairs are able to share temporal coherence over the entire chip, 
which would significantly improve the quality of quantum functionalities built upon. 

The superior performance of the demonstrated device offers a great opportunity for realizing 
a variety of high-quality quantum functionalities on chip. Straightforward applications include 
producing heralded single photond^ and quantum entanglement in various degrees of freedom 
(frequenc)^^, pathP, timeP, polarizatiorl^, etc.), which are essential for quantum communication 
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and quantum computin^^Sl Compared with conventional approaches, the high photon coherence 
in the device would significantly improve the purity of the constructed quantum states. As the 
cavity enhancement only occurs within the cavity photon lifetime, the demonstrated approach is 
able to operate with a repetition rate above gigahertz by use of pulsed pumping. Further increase 
of photon-pair coherence and spectral brightness can be obtained by improving optical Q (e.g., 
through surface passivatiorP'). Such a high coherence would enable teleporting quantum entan- 
glement intra-chip, inter-chip, or off-chip over far distance free from the stringent synchronization 
requirement in conventional scheme^^J^Mil ^jj^g direct photon coupling into optical fibers with 
negligible losses (Fig.|2]) is particularly suitable for this application). Moreover, the comb nature 
of emitted photons makes readily available wavelength-division multiplexing in quantum regime, 
and thus significantly extends the spectral efficiency and wavelength management capability of 
integrated quantum photonics. In particular, it exhibits great potential for realizing multipartite 
entanglement^ and cluster state"^ inside a single device, which would dramatically enhance the 
capability of integrated quantum computing. On the other hand, the high coherence and spectral 
flexibility of the demonstrated approach enable excellent interfacing with other chip-scale infor- 
mation processing/storage element^. These features together with the developed silicon photonic 
technolog)!™ 

offer a great opportunity to ultimately form a complete CMOS compatible chip- 
scale platform for quantum information processing. 

Methods 

The silicon microdisk resonator is fabricated from a standard silicon-on-insulator wafer by use of e-beam lithogra- 
phy to define the device pattern, fluorine-based plasma to etch the silicon layer, and hydrofluoric acid to undercut the 
silica pedestal. Light is coupled into and out of the microdisk through a tapered optical fiber which is anchored onto 
two nanoforks fabricated near the microdisk for stable operation (see Fig. |2]i. Touching to the nanoforks introduces 
a small insertion loss of only about 7%, as indicated by the cavity transmission spectrum given in Fig. |3^. The PL 
spectra were recorded at the transmission ports of the CWDM DEMUX to suppress the residual pump wave which 
otherwise would cause significant stray light inside the spectrometer. 

The CAR is obtained by the following procedure: First, at each pump power, a delay-time-dependent coincidence 
counting histogram (we term as coincidence spectrum here for simplicity. See Fig.[5]for an example) is obtained with 
zero time delay between the two SPDs. The coincidence spectrum is integrated over its full width at half maximum 
T^FWHM to obtain the total coincidence C. The statistical error of C is given by Ec = VC since the photon counting 
follows a Poisson distribution. To find the accidental coincidence, a series of coincidence spectra are obtained with 
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a time delay of JTq between the two SPDs, where Tq is the detector clock period and /' = 1,2, ...,200. Integrate each 
of them within the same TfWHM to obtain an accidental coincidence Aj. The average and the standard deviation of Aj 
give the accidental coincidence A and its statistical error Ea- Note that no any accidental coincidence is subtracted. 
Finally, the CAR is given by CAR = (C - A)/A and its statistical error ^car is given by \J (^^^ + (j^^ . 

To reduce the statistical error, each data point in Fig. |4] was recorded with a time period of 5400 — 18900 seconds 
depending on the magnitude of photon-pair flux. With such amount of photon counting data, ^ is very small (< 2%) 
and is dominantly determined by ^ since the accidental coincidence is very small in our device. 

To find the photon-pair emission flux shown in Fig.|4^-h, a same procedure is carried out but now the integration is 
performed over the entire time window of coincidence spectra (20 ns, with a detector gating width of 10 ns) to obtain 
the total coincidence Q, the accidental coincidence At, and their statistical errors Eq and Ea,- The true coincidence 
is thus given by Ct — Cr — At. Its statistical error Ecj^ is given by Ecj = \J (Eq)^ + (Ea,)^- Dividing Ct by the total 
number of detector gates within the data collection time we find the true coincidence counts per gate (which is shown 
on the right axes of Fig. [4^-h). The photon-pair emission flux is obtained by calibrating the true coincidence counts 
per gate with the clock frequency, duty cycle, and quantum efficiency of the detectors as well as the propagation losses 
from the device to the detectors for the photon pairs. The statistical error of pair emission flux is obtained accordingly 
from Ecj^ . 

The IRF of our coincidence counting system is obtained with identical ultrashort pulses launched onto the two 
SPDs. The ultrashort pulses come from a mode-locked fiber laser with a pulse width (FWHM) of 240 fs. They are 
broadened to about 5 ps by the delivery fibers before they hit on the SPDs. The pulse width is much shorter than the 
timing jitters of the SPDs, thus providing a fairly accurate approach to characterize the IRF. 

In all photon counting measurements, the detector gating width is set to 10 ns and the deadtime is equal to the 
clock period. The two InGaAs SPDs have dark counts of 5.09 x 10^^ and 3.18 x 10^^ per gate, respectively, at a 
quantum efficiency of 15% and a clock frequency of 250 kHz. The detector dark counts increase to 8.96 x 10^^ and 
6.93 X 10^^ per gate at a quantum efficiency of 25% and a clock frequency of 125 kHz. In both cases, the after-pulsing 
probability is less than 8%. 
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Appendix A: Group-velocity dispersion of a silicon microdisk resonator 

For each mode family, the silicon midrodisk resonator is simulated by the finite-element method 
to find the cavity resonances (0(m) as a function of mode number m over a broad spectrum. As 
the propagation constant of a cavity mode is well approximated hy k = m/R where R is the device 
radius, we obtain the dispersion relation k{o)) for the cavity modes. Fitting it with a high-order 
polynomial, we obtain the group-velocity dispersion shown in Fig. lb of the main text. 

In general, the group-velocity dispersion of a microdisk resonator is close to a slab waveguide 
with the same thickness, which can be used as a rough guidance for searching the group-velocity 
dispersion of a microdisk resonator. 

Appendix B: Theory of photon-pair generation in a microdisk resonator 

In this section, we provide the theory describing photon pair generation in a silicon microdisk 
resonator. We consider only the Kerr nonlinear optical interaction and neglect two-photon absorp- 
tion and free-carrier effect because of their negligible effect in our device. We consider first the 
case shown in Fig. |6] where a pump wave launched into the resonator produces a pair of signal and 
idler photon which are transmitted out into the coupling waveguide. The Kerr nonlinear interaction 
among the three cavity modes at frequencies CQoj U = Pi^J) can be described by a Hamiltonian 
H = Hq+Hj, where Hq describes the passive cavity modes coupled to the bus waveguide and Hj 
describes the Kerr nonlinear interactiorP^ inside the cavity, with the following forms 



Ho= {fiO}oja]aj-h./r^j a]bje-'"'j' +b]aje''^j' }, (Bl) 

j=P,s,i 

Hi = -^gp{al)^al - lhalap (^gpsalas + gpiaja,^ - hgpsialajal - hgl^^{alfasai, (B2) 



where the intracavity field operator aj (j = p,sj) is normalized such that a jay represents the 
photon number operator and bj is the field operator of the incoming wave at carrier frequency 
COj inside the coupling waveguide normalized such that b^jbj represents the operator of the input 
photon flux. {b%) = {b]bi) = since only the pump is launched into the cavity, bj satisfies the 
commutation relation of [bj{t),b^j{t')] = d{t — t'). Egj is the external coupling rate of the cavity 
mode at (Ooj. 



In Eq. (B2), the first and second terms describe the self-phase and cross-phase modulation 
(SPM and XPM) from the pump mode, respectively, and the third and fourth terms govern the 
four- wave mixing (FWM) process. We have assumed the signal and idler are much weaker than 
the pump so that we can neglect the SPM and XPM introduced by them, gp, gp^ (j = i), and gp^i 
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FIG. 6: Schematic of photon pair generation in a silicon microdisk resonator. 

are the vacuum couphng rate for the SPM, XPM, and FWM, respectively. They are given by the 
following general expression 



,?ijkl 



3^ ^ ftlOi £L»Oj £t»Ok ft>01 T?ijkl (3) , , 



(B3) 



4eon,nj%«z^jki 

where n,, (v = i,j,k,l) is the refractive index of silicon material at COov and X^^^ is the third-order 
nonlinear susceptibility, t^jki = {ViVjVkViy^^ is the average effective mode volume and Vy (v = 
i, 7, k, I) is that at individual frequency COov given as!^ 

2 



Jdrer{r,(Ooy)\Ey{r,(Oo^)\^^ 



(B4) 



J^idr£^{r,cOo^)\Ey{r,cOoy)\'^ 

where er{r, COov) is the dielectric constant of material and the subscript si in the integral denotes 
the integral over the silicon core. T/ijki describes the spatial mode overlap given 



Y{v^i,i.k,i!sidr£}y\Ey\ 



1/4 



(B5) 



where Crv = £r{f: COov) and Ey = Ey(r, (Oov) is the electrical field profile for the cavity resonant 
mode at (Oov- 



The vacuum coupling rates for various nonlinear effects are thus obtained from Eq. (B3) as 
gp = ,?pppp, ,?pj = ^pjjp U = i), and ^psi = gspip- In our device, the three cavity modes are of close 
frequencies and similar mode profiles, resulting in 

cr\n2h(0py/0)s(0i 



gp ~ gpj ~ ^psi = g 



TlsTliV 



(B6) 
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where 7] and V denote Tj^pip and Vspip^ respectively, and we have used the definition of Kerr non- 



at the pump frequency. Using Eqs. (Bl I and (B2 1 and counting in 



hnear coefficient^^: n2 

the intrinsic cavity losP, we obtain the following equations of motion in the Heisenberg picture 
governing the wave dynamics inside the cavity: 



dt 
das 

dt 
dui 

dt 



i(Oop - rtp/2)ap + iga a + i^yTepbpe ' p +i^JTQpUp, 



(B7) 



{-i(OQ^-Yl^/2)as + 2iga^papas + iga]ap + i^JY^J)se +i^/rosUs, (B8) 
{-iCOoi - Tii/2)ai + 2igalapai + igala^p + i^/Y\ibie^'^'' + z a/T^m,-, (B9) 



where Foj and = Foj +rej (j = p,s,i) are the photon decay rates of the intrinsic and loaded 
cavity, respectively. My is the noise operator associated with intrinsic cavity loss, which satisfies 



the commutation relation of [uj{t),u^j{t')] = 5{t — t'). In Eqs. (B7)-(B9), we have neglected the 
pump depletion and the signal/idler-induced SPM and XPM. The transmitted field is given bjJ^ 



fj = bj + iy^jaj. 



(BIO) 



In general, the pump mode can be treated as a classical field and Eqs. (B7)-(B9) can be easily 



solved to find the solutions for and a,, from which, together with Eq. (BIO), we find the emitted 
photon fluxes of signal and idler as 

Tesrti 



iflmit)) 



2% 

2% J- 



\Bi(o)\ d(0. 



\B(a})\ dco. 



and the pair correlation as 



PcitsA) ^ {fKt.)fKts)fs{ts)md) - {fi{ts)Mts)){fjm{td) 



— r r 



2n 



B{-(o) [rtsA((o) - 1] e-"'^'-^-"U(0 

-oo 



where A(tt)) and 5(a)) have the following expressions 

rti/2-zCi) 



A (0) 



B{co) 



|2' 



{rj2 -ico){r,,/2~ico)-{gNp 

-iS^l 

{r,,/2- ico){r,,/2-ico)-{gNp)^-' 



(Bll) 
(B12) 



(B13) 



(B14) 
(B15) 



where Np = {a^ap) is the average photon number of the pump wave inside the cavity. 
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Equations (|B1 1|)-([BT3]) are quite general since they include the multiphoton generation induced 



by the stimulated FWM. For single photon-pair generation, Eqs. (Bl 1 ) and (B12) reduce to 



{fKt)fs{t)) 

iflmit)) 



(B16) 
(B17) 



where T = (Fts + rti)/2 represents the average photon decay rate of the loaded cavity. The photon- 
pair correlation reduces to 



r2 



f2 



KgNp 



"P) 



if ts > ti 
if ts < ti 



(B18) 



Equation (B18) can be regarded as the joint probability density of emitting a photon pair into the 
coupling waveguide. It is given as Eq. (1) in the main text. Integrating Eq. (B18) over T = ts — ti, 
we obtain the average photon-pair emission flux as 



Rr 



FesFei 2{gNp 



(B19) 



TtsFt, F 

which is Eq. (2) in the main text. Comparing Eq. ( B19[ ) with Eqs. (B16) and (B17), we can 
find that the photon-pair generation rate inside the cavity is 2{gNp)^/T and the photon extraction 
efficiency is Fej/Ftj (j = sj) for the individual signal and idler mode while it is FesFgi/ (FtsFti) 
for the correlated photon pair. As Fgj/Ftj is less than 1, the photon-pair flux is smaller than the 
individual photon fluxes. This is experimentally observed in Fig. 4 in the main text. 

In practice, the situation is much more complicated than the ideal situation described above, 
since, due to the light scattering, the optical wave inside the cavity will couple to the degenerate 
mode propagating backward in the opposite direction (Fig. |7]). For a passive cavity, this phe- 
nomenon manifests as a doublet in the transmission spectrum (Insets of Fig. 3a in the main text). 
In this case. The Hamiltonian Hq and Hj become 



Hi= («Jf)^«pf+(«Jb)^«pb+4«Jf«pf^^pb«pb 



j=p,s,i 

-njr 



(B20) 



psi 



-t ^2 



t ^2 



'^sf^if'^pf + '^sb'^ib'^pb 



- %psi («Jf)^«sf«if + («pb)^«sb«ib 



t \2, 



(B21) 
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FIG. 7: Schematic of photon pair generation in a silicon microdisk resonator. The clockwise and counter-clockwise 
modes are coupled inside the cavity. 



where ajf and ajb (j = p^s^i) represent the optical fields propagating clockwise and counter- 
clockwise, respectively, inside the cavity (Fig.[7]), and j8j is the coupling coefficient between them. 
bjf and Z?jb represent the input fields propagating forward (input from the left end) and backward 
(input from the right end), respectively, inside the coupling waveguide. The transmitted fields now 



become fj = b^f + i^y^Qja^f. Using Eqs. (B20) and (B21 1 and following the same procedure, we 
can find the equations of motion as 



da 



pf 



dt 
da^b 

dt 
da^f 

dt 
dttsh 

dt 
dan 



(-/cODp - rtp/2)apf + i^pQ^h + ''^(«pf«pf + 2aJ(,apb)apf + /^pf , (B22) 

(-/(OOp - rtp/2)apb + //3papf + ig{a^^^a^y, + 2ajfapf)apb + /^pb, (B23) 

(-zCODs - rts/2)asf + //3.,asb + 2ig{a^^^a^f + a^^^a^h)aif + iga-^a^f + z'Csf , (B24) 

(-zcODs - rts/2)asb + /Ar«sf + 2z^(ajfapf + aj^apb)asb + «^«ib^pb + 'Csb, (B25) 



dt 
dci[\) 

dt 



(B26) 
(B27) 



'-i(Ooi - rti/2)aif + z'Aaib + 2zg(ajfapf + aj^apb)aif + igal^a^f + z'Cif , 

[-iOoi - rti/2)aib + z/3*aif + 2ig{a^p^apf + flp^,apb)aib + «^«sb^?b + ^'Cib 

where Cjv = y^i^jye^'^^' + A/Tq^^jv (j = i and v = f,b). 

Similar to the previous case, Eqs. (B22)-( [B27] ) can be used to find the solution of signal and 
idler fields inside the cavity ajv (j = s, i and v — f, b), through which we obtain the signal and idler 
fluxes as 

{fhmt)) 



iflmit)) 



resrti 

In 
In 



|Mi3(G))P+|Mi4((0)P 



|M3i((0)p + |M32(a))|' 



dco, 
dco, 



(B28) 
(B29) 
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and the pair correlation as 



Pc(t„ti) = ifH'dflMfsMM'i)) - {flMf,M){fI{'dM'i)) 

In J-o 



(B30) 



where Mj^ is the elements of the matrix M{(o) which is given by 











-//3, /ft)-rti/2/ 



(B31) 



Equations (B28) and ( B29| ) provide the theoretical curves of individual signal and idler fluxes 



given in Fig. 4a-d in the main text. Equation (B30) shows that the photon-pair correlation is only 
a function of time difference T = f. , — PcitsiU) = Pc{t^)- Integrated pd'c) over T, we obtain the 
photon-pair emission flux Rc = J^^ Pc{^)dT, which is given as the theoretical curves in Fig. 4e-h 



in the main text. The complete expressions of Eq. (B3)-(B5) are used in the theoretical calcula- 
tion, where the effective mode volumes and the mode overlap factors are obtained from the finite 
element simulation of the device. A Kerr nonlinear coefficient, n2 = 5.0 x 10^^ cm^/GW, is used 
in the theoretical calculation. The cavity parameters, such as Foj, Fej, and /3j (j = p,sj), are 
measured from the laser-scanned transmission spectrum of the passive cavity. The Pi mode at 
1584.1 nm is outside the scanning range of our lasers and its cavity parameters cannot be obtained 
through independent measurement. We approximate its intrinsic Q and mode splitting with those 
of the nearest Pi mode at 1564.9 nm since the properties of cavity modes remain similar within 
each mode family. Moreover, we approximate its external coupling with that of the P2 mode at 
1553.6 nm because of their similar external couplings as observed from the PL spectra shown in 
Fig. 3b and c of the main text. Figure 4d,h,i of the main text show that these parameters provide a 
good approximation for the cavity properties of the Pi mode at 1584.1 nm. 



Appendix C: Statistics of coincidence photon counting 

In this section, we provide the theoretical description of the statistics of coincidence photon 
counting with gated single photon detectors (SPDs). Assume a signal and idler photon strike the 
detectors at time ts and ti with a probability density of PphitsJi)- The SPDs then produce a pair 
of photoelectric events which are recorded by the coincidence counter at a later time ti = ts + Ti 
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and ?2 = ti + "^2 (Fig- [8]), where Ti and T2 include the electron transit time inside the SPDs, the 
propagation time from the SPDs to the coincidence counter, and time delays introduced by the 
coincidence counter, with a joint probability density of py(Ti, 12)- The joint probability density of 
the photoelectric events is thus given by 

<■ /.+00 

PSl{ti,t2) = ^1^2 // Pph{ts,ti)pj{X\,X2)dXidX2 

J J — 00 

Pph{t2-ti - TD)pmp{XD)dTD, (CI) 

-00 



where t,j (j =1,2) are the detection efficiencies of the two SPDs, and we have used the 
fact that Pph{tsJi) = Pph{ti — h) which is only a function of the time difference between the 



signal and idler arrival. For a correlated photon pair, Eq. (B30) indicates that Pph{tsJi) = 
asOCi{fJ{ti)f^{ts)fs{ts)fi{ti)) where as and a,- are the propagation losses of signal and idler pho- 



tons from the coupling waveguide (a tapered fiber in our case) to the SPDs. In Eq. (CI ), /'irf(t) 
is the instrument response function (IRE) of the coincidence counting system, given as 

/'IRf(Td)=/ py(Ti,Ti+Ti))JTi. (C2) 

Jo 

The IRE is independently measured by launching identical ultrashort pulses onto the SPDs and 
recording the coincidence spectrum (see Method in the main text). In out coincidence counting 
system, the IRE is dominantly induced by the timing jitters of the two SPDs. As the detector timing 
jitter varies with its quantum efficiency (QE), the IRE is recorded at different QEs for specific 
experiments (QE is 15% and 25% for Eig. 4 and 5 in the main text, respectively). In general, the 
experimentally recorded IRE is well described by a Gaussian function since the timing jitters of 
the SPDs follow a Gaussian distribution. 

As a detection event can be either a photoelectric event fired by a photon or a detector dark 
count, the joint probability density, p{ti , ^2), that the two detectors produce a pair of events at time 
ti and t2 (Eig. [8]) is thus given by 

p{hj2) = Psiih.h) + Psoituh) + PDi{hJ2) + PDD{tut2), (C3) 

where pso, Pdi, and poo are the joint probability densities of an event pair of signal— dark count, 
dark count— idler, and dark count— dark count, respectively, produced in the two separate SPDs. 
The generation of a detector dark count is independent of an event generated in another detector 
(either a photoelectric event produced by an incoming photon or a detector dark count). As a 

result, PSD{hj2) = Ps{tl)PD2{t2), PDlitlJl) = PD\{h)Pl{t2), and PDD{tht2) = PDl{tl)PD2{t2), 
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FIG. 8: a. Schematic of coincidence photon counting. A signal and idler photon arrive at the detectors at time and 
f,, respectively, with a probability density of Pph{ts,ti)- The two detectors produce a pair of events which are recorded 
by the coincidence counter at time fi and t2, with a probability density of p{t\,t2)- b. The time relationship of the 
events arriving at the coincidence counter Each event arrives within a gating window centered at tjQ (j ~ 1, 2) and 
with a width of 7g. The two input channels of the coincidence counter has a time delay of 7j. 

where ps{t), Pi{t), and pojit) U = 1,2) are the probability densities of detection events from 
signal, idler, and dark counts, respectively, produced in the two detectors. In experiment, ps{t), 
pi{t), and PDjit) are equal to the average detection rates of the related events, which can be 
measured independently. Consequently, Eq. (|C3[) becomes 



p{tuh)= Psiih Psih ) Pd2 (tl) + Pdi ih ) Pi (ti) + Pdi (h ) Pdi (h)- (C4) 

Therefore, the coincidence dXx = t2 — t\ within a small delay-time interval of dx is given by 
p{z)dx= / dti / dt2p{ti,t2) 

Jtio-Tc/2 Jt2o-Tc/2 
= (To - I T - Trf I) [j35/(t - Td) + PSPD2 + PDlPl + PD\PD2] , (C5) 



for \x — Td\ < Tq and p(t) =0 otherwise. In experiment, dx corresponds to the time resolution 
of the coincidence counter, which is 4 ps in our system, tio and ?20 the centers of the gating 
windows in the two channels (Fig. [8]). As a result, the delay-time-dependent coincidence counting 
histogram (we term as the coincidence spectrum here for simplicity) is obtained accordingly as 



C{x)dx = RGTp{x)dx, 



(C6) 



where Rg is the clock frequency of the SPDs and T is the data collection time. 



Equations (jCTj), ( |C5| ), and ( |C6[ ) provide a complete description of coincidence photon counting, 
from which we can obtain the true coincidence spectrum Ct{x) and the accidence coincidence 
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spectrum Ca(t) as 



Cr(T) Jt = RcTdx (Tg - |t - r^l) [psi{x - T^) - psPi] 

/+°° 
Pc{t -Td- TD)piRFMdrD, (C7) 
-oo 

Ca {x)dx = RcTdx iTG-\x-Td\){ps + PDi)iPi + PD2), (C8) 



where pdx) is given in Eq. (B30). The theoretical coincidence-to-accidental ratio (CAR) shown 



as the red curves in Fig. 4e-h of the main text is obtained as 

-7fwhm/2 . 



CAR = ^_ ^^^^ 

f-TZ:l2CAi^ + WdX 



where Tfwhm is the full width at half maximum of the coincidence spectrum. In Eqs. (C7) 

and H), PS = +m and pi = ^20ti{fJ{t)fi{t)) + pr^, where and 



{fi{t)fi{t)) are given by Eqs. (B28) and (B29), and pR^ and pr^ are the probability densities of 



detection events from Raman noise photons at the signal and idler channels, respectively (See 



the next section for the details of characterizing Raman noises). On the other hand, Eq. (C8) 
shows clearly that the accidental coincidence can be accurately calculated from independently 
measured ps, pi, and poj (j = 1 , 2). Therefore, using these measured parameters together with in- 
dependently recorded true coincidence we can reproduce the experimentally recorded CAR. This 
theoretical approach explains well the saturation of CAR at low pair fluxes and is commonly used 
in literature!^. It is shown as the black dashed curves in Fig. 4e-h of the main text. 

The theoretical curves shown in Fig. 5 of the main text are described by Eq. ( |C6| ) (normalized 
by its peak value), in which Pphits.ti) = (Xs(Xi{f^{ti)f^{ts)fs{ts)fi{ti)) obtained from Eq. (B30) and 



j9irf(t) is given by a Gaussian fitting to the experimentally recorded IRE. The cavity parameters 
used in {fj {ti)fj {ts)fs{ts)fi{ti)), such as Foj, Fej, and (j = p, s, i), were independently measured 
from the laser-scanned transmission spectrum of the passive cavity. 

Appendix D: Characterization of Raman noise produced in delivery fibers 

To characterize the noise for each photon mode, we recorded the photon flux when the center 
wavelength of the tunable bandpass filter was tuned away from the photon mode so that the latter 
was completely blocked (see Fig.|9]). The pumping condition remained the same as what we used 
to characterize the photon pairs, with the tapered fiber touched down to the nanoforks and the 
pump wavelength on the cavity resonance. Such an experimental condition ensures that the pump 



22 



Pump 



fiber 



'r< 


~5m 




~5m j 






— H< — 


>; 



fiber 1 ^^^^ fiber 2 




signal 



fiber 3 J 
fiber 3 




CWDM CWDM '^'^"^ 

MUX Si microdisk DEMUX „ r ,. 

Bandpass Filter 



SPD1 
SPD2 



FIG. 9: Schematic of the experimental setup. Details can be found in Fig. 2 of the main text. 

power distribution inside the delivery fibers remains unchanged. The photon flux recorded under 
this condition is not produced from the device since the photon wavelength does not coincide with 
any cavity resonance. 



Figure 10 shows an example of recorded photon fluxes for Pj, where the solid and open circles 



show the cases when the filter centers were tuned onto and away from the photon mode pair, 



respectively. The noise photon fluxes exhibit a linear dependence on the pump power. Figure 10 
shows clearly that the observed noise is the dominant noise source accompanied with the photon 
pairs. To further verify that the observed noise does not come from the device, we directly lifted 
the taper up so that the tapered fiber did not couple to the device and the pump wave propagated 
only inside the delivery fibers (as well as the optical components in the experimental system. See 
Fig. |9]). Kept other conditions unchanged and recorded again the photon fluxes, we obtained the 



open triangles shown in Fig. 10 which are very similar to the open circles. This clearly verifies 
that the noise photons are not generated from the device but produced inside the experimental 
testing system. 

In our testing system (Fig. |9]), the CWDM MUX has a band isolation greater than 120 dB. 
As a result, any broadband noise produced before the MUX (e.g., in fiber 0) will be completely 
blocked. On the other hand, the CWDM DEMUX separates the pump wave from the photon 
pairs and no noise will be produced m fiber 3 after the DEMUX. Moreover, when the pump wave 
is nearly critically coupled into the device, the pump power in fiber 2 after the device is small. 
Therefore, the noise should be dominantly produced inside fiber 1 before the device, which has a 



rough length of ~5 m. The open triangles in Fig. 10 have higher magnitudes than the open circles 



simply because, when the taper is lifted up, the pump power will increase in fiber 2 and will start 
to produce noise. 



Figure 10 shows that the idler channel at the Stokes side of the pump has a larger noise photon 



flux than the signal channel at the anti-Stokes side. This feature together with their linear power 
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FIG. 10: Photon fluxes recorded for Pj, wifli detector dark counts subtracted. The solid and open circles show the 
photon fluxes when the center wavelengths of the bandpass filters are tuned onto and away from the Pj mode pair, 
respectively, with the tapered fiber touched down to the nanoforks. The open triangles are recorded when the tapered 
fiber is lifted up and thus does not couple to the device. The blue and red show the signal at shorter wavelength and 
idler at longer wavelength, respectively. The dashed lines are only for eye guidance. 

dependence implies that the noise is likely to be the Raman noise produced inside fiber 1 . To 
verify this, we incremented the length oi fiber 1 (with the tapered fiber not coupled to the device) 
and monitored the noise spectrum. Figure [TT^ shows an example of the noise spectra recorded 
with three different fiber lengths of ~5 m, ~29 m, and ~1 km. The similarity of the spectra 
among these three cases indicates clearly that they share an exactly same physical origin. The 
noise spectra show the Raman noise characteristics of silicaF', with the Boson peak clearly visible 
at a frequency shift of ~ 1.5 THz away from the pump. 

To reveal more characteristics of the noise spectra, we tuned the pump wavelength to the left 
end of the CWDM and recorded again the noise spectra (Fig. [TTJ5). In this case, the noise spectra 
primarily show the Stokes side. Both the Boson peak and Raman gain peak are clearly visible, with 
a frequency shift of about 1.5 and 13 THz, respectively, away from the pump frequency. Tuning 
the pump wavelength to 1570nm we recorded the noise spectra now mostly on the anti-Stokes 
side (Fig.pT|:). Again, we observed both peaks but with smaller amplitudes. These observations 
clearly verify that the observed noise indeed comes from the spontaneous Raman scattering inside 
the delivery fiber. 
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As the Raman noise is very broadband, in experiment, we used a pair of tunable bandpass 
filters with a 3-dB bandwidth of 1 .2 nm in front of the SPDs (Fig. |9]) to cut the Raman noise. The 
Raman noise fluxes associated with each mode pair were measured with the filter centers tuned 
away from the mode pair far enough that the latter were completely blocked. The recorded Raman 



noise fluxes were used in Eq. (C8) for finding the theoretical CAR, as discussed in the previous 
section. Due to the strong frequency dependence of Raman noise spectrum, a certain uncertainty 
may be accompanied with the recorded Raman noise fluxes. Note that the Raman noise can be 
further suppressed in the future by using narrower filters since the correlated photon pairs have 
much narrower spectra. The noise can also be reduced by shortening the length of fiber 1. Such 
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FIG. 11: a. Raman noise spectra recorded with different lengths of fiber 1, with an input pump power of about 
270 IJ.W. The spectra were recorded at different transmission ports of the CWDM DEMUX (indicated by color), for 
easy suppression of the pump. The solid curves and dotted curves show the cases with a fiber length of ^29 m and 
~1 km, respectively. The amplitudes of the dotted curves were reduced by a factor of 10, for easy comparison with 
the thin solid curves. The inset shows the case with a fiber length of ^5 m, same as what we used to characterize 
correlated photon pairs. The pump power and wavelength are the same for these three cases. The spectral cut-off 
beyond ~1590nm is due to the quantum-efficiency cut-off of the detector used in our sepctrometer b. Raman noise 
spectra with the pump located at 1470nm, showing the Stokes side, fiber 1 has a length of ~' 1 km. c. Raman noise 
spectra with the pump located at 1570nm, showing mostly the anti-Stokes side. Fiber 1 has a length of ^1 km. 
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Raman noise will be completely absent if an on-chip lasei is used as the pump. 
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